Stochastic Acceleration of Low Energy Electrons in Cold 

Plasmas. 



Julia M. Pryadko 1 and Vahe Petrosian 2 

Center for Space Science and Astrophysics 

Stanford University 

)£ ■ Stanford, CA 94305-4055 

ON 
ON 

ABSTRACT 

O 

O 

We investigate the possibility of stochastic acceleration of background 

low-energy electrons by turbulent plasma waves. We consider the resonant 

interaction of the charged particles with all branches of the transverse plasma 

waves propagating parallel to a uniform magnetic field. Numerical results and 

asymptotic analytic solutions valid at non-relativistic and ultra-relativistic 

energies are obtained for the acceleration and scattering times of electrons. 

These times have a strong dependence on plasma parameter a = u pe /Q e (the 

ratio of electron plasma frequency to electron gyrofrequency) and on the spectral 

^ • index of plasma turbulence. It is shown that particles with energies above 

Q-T 

certain critical value may interact with higher frequency electromagnetic plasma 
waves and this interaction is allowed only in plasmas with a < 1. We show 
that for non-relativistic and semi-relativistic electrons in low-a plasmas the 
ratio of the acceleration time to the scattering time can be less than unity for 
a wide range of energies. From this we conclude that the transport equation 
derived for cosmic rays which requires this ratio to be much larger than one 
is not applicable at these energies. An approximate "critical" value of particle 
energy above which the dynamics of charged particles may be described by 
this transport equation is determined as a function of plasma parameters. We 
propose new transport equation for the opposite limit (energies less than this 
critical value) when the acceleration rate is much faster than the pitch angle 
scattering rate. This equation is needed to describe the electron dynamics in 
plasmas with a <; 0.1. 
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1. INTRODUCTION 

The acceleration and propagation of charged particles in magnetized plasmas via their 
stochastic interaction with turbulent plasma waves is a problem of wide interest in different 
astrophysical areas. The interaction of the particles with plasma turbulence causes diffusion 
of the energetic particles in phase space and leads to stochastic acceleration which is a 
second order Fermi process. The rates of these interactions are controlled by the magnetic 
field B Q and other physical properties of the background plasma such as the background 
particle density n, the energy density and spectrum of the turbulence and, if in thermal 
equilibrium, the background temperature T. 

The influence of different modes of plasma waves on the dynamics of high energy 
charged particles has been discussed in the literature for different plasma conditions. 
Whistler waves for electrons in an electron-proton background plasma for solar flare 
conditions have been considered by Steinacker and Miller (1992), Miller (1991), Miller 
and Ramaty (1987), Hamilton and Petrosian (1992) and by Achatz et al (1993) for 
interplanetary conditions. Benz (1977) and Melrose (1980) consider Langmuir (plasma) 
waves, Achterberg (1981) and Beaujardiere and Zweiles (1989) deal with magnetosonic 
waves, Benz and Smith (1987) describe interaction with lower hybrid waves. The transport 
equation for transverse Alfven waves in a cold background medium has been considered by 
Schlickeiser (1989) and this consideration have been extended for electromagnetic branch 
in general astrophysical plasmas by Dung and Petrosian (1994) (hereafter refered to as 
DP). Most of these and other works investigating particle acceleration deal with relativistic 
particles and assume an initial source of non-thermal high energy particles. 

We are interested in the investigation of acceleration of the particles present in the 
background. These particles usually have a Maxwellian distribution with kT <C mc 2 and 
therefore have low energies. The general problem is formulated in DP and some results 
are given for high and intermediate energies. In this paper we continue this investigation 
of interaction with all modes of transverse plasma waves propagating parallel and/or 
antiparallel to the ordered magnetic field with emphasize on low energy electrons. In §2 
we give an overview of the basic equations of turbulent plasma theory used in this paper 
and define a ratio of scattering to acceleration times in terms of Fokker-Planck coefficients 
which are discussed in §3. In §4 we evaluate this ratio for different plasmas and particles of 
different energies and suggest a new transport scenario for low energy particles. In §5 and 
§6 we derive approximate analytic expressions for the above ratio and the acceleration time 
for non-relativistic and extremely relativistic electrons, respectively. A brief summary and 
our conclusions are presented in §7. 
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2. BASIC EQUATIONS 

We consider the behavior of the energetic charged particles in a background plasma 
with uniform magnetic field of strength B and superposed plasma waves. The gyrophase- 
averaged phase space density f(z,t,p,fi) then obeys the differential Fokker-Planck equation 
(Schlickeiser 1989): 

df df 1 d f d d\ d ( d d\ 

m +Vfi d- z = J 2 ^/ [ Dpp d-p + Dp "dH ) f + dH + D ^d~ P ) f + s - (1) 

Here z is the distance along the field lines, S is the source function, /x denotes the 
cosine of the pitch angle, v = p/irwy) is the velocity and 7 = (1 + p 2 / [mc] 2 ) 1 ^ 2 is the 
Lorentz factor of the particle with momentum p. The interaction of the particles with 
plasma waves is described by the three Fokker-Plank coefficients D^ p = D Ptl and D pp 
which depend on the properties of the turbulence. 

Solving the differential equation ([[]) in general is not simple so one usually considers 
the solution for certain physical conditions when it can be simplified. In cases when 
the pitch angle scattering time scale T sc (/i) ~ £)~* is much shorter than the traverse 
time Ttr — L/v, where L is the size of the turbulent plasma region, the pitch angle 
distribution will be nearly isotropic so that the anisotropic part of the phase space density 
g(z,t,p,pi) = f(z,t,p,/j l ) — F(z,t,p) is much smaller than the pitch angle averaged phase 
space density 

F(z,t,p) = \ ( + dfif{z,t,p,fi). (2) 



2 J- 
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In addition if the ratios 



D pp/P 2 // 1 T) 1 ,, „\ D p^P 



RiM = -TT~ « x > **M = -jt^ « 1 (3) 

MM MM 

then equation (|l|) can be reduced to the diffusion-convection equation which is also known 
as the transport equation for the cosmic rays (Jokipii 1966; Kirk et al. 1988; Schlickeiser 
1989; DP): 

— - — K — + Uyv) — — - — — (p 3 vk )— + — — (A — + Q(z t p) (4) 

dt dz dz dz dp p 2 dp dz p 2 dp dp 

Here Q(z,t,p) is the pitch-angle averaged source term and the three transport coefficients 
are defined as 

v 2 f +1 , (l~/i 2 ) 2 
8 J-i 
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1 r+ L 

K2 = 4 J i _ ^ ^ 2 ^ 

K3 = l£dLiD^(R 1 -R 2 2 ) (7) 

The first four terms of the right hand side of equation (ffl) represent spatial diffusion, 
adiabatic acceleration/deceleration, spatial convection and momentum diffusion. If the 
above equations are satisfied we can define averaged scattering, acceleration and spatial 
diffusion times as t sc = 8ki/v 2 , r ac = 1/k, 3 , r di ff = L 2 / K\ = 8T 2 r /r sc , respectively. 

The assumptions that lead to transport equation (^) are not always valid. While the 
condition t sc <C L/v holds for wide range of particle energy in most astrophysical plasmas 
the requirements in equation (|3|) are not always satisfied. As shown below these ratios can 
exceed unity for low energy background particles. 

The goal of this paper is to investigate the range of validity of equation @ and suggest 
a different transport scenario for cases when this equation is inapplicable. Here we will deal 
only with waves propagating parallel or antiparallel to the magnetic field. The influence of 
transverse plasma waves propagating perpendicular to the magnetic field on the dynamics 
of charged particles will be considered in a subsequent paper. 



3. FOKKER-PLANCK COEFFICIENTS 

Following DP (see also Lerche 1968) and assuming a power law distribution of plasma 
turbulent energy density as a function of wave vector, S(k) — (q — l)S to t Kmin K~ 9 (f° r 
K > K m i n and q > 1), it can be shown that the Fokker-Planck coefficients necessary for 
evaluation of the ratios in equation (0) can be written as: 

D m = ^(1 - £ (l - ^^f 1 ) 2 X(k 3 ) (8) 

3* = J- n _ ^ £ MM L _ x[k ) (9) 
p t pi pi P V P J 



where 



p2 T ^i-tf£^Uft). do) 



\K 3 \ 
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(3 = v/c and the dimensionless wave vector kj is one of the roots (maximum of four) of the 
resonant condition: 

u(kj) - n/3kj ^1/7 = 0, kj = Kjc/n e . (12) 

Here and in what follows the upper and lower signs refer to the right (R) and left(L) hand 
polarized plasma modes. The wave frequency Q, in units of electron gyrofrequency fl e , is 
determined from the dispersion relation: 

£ =i - ( *S(L±f) - ^ )=n « (13) 

where 5 = — is the ratio of electron to proton masses and 

a = u pe /VL e = 3.2 K/10 10 cm" 3 ) 1 / 2 (5/ 100 G)" 1 (14) 

is the ratio of electron plasma frequency to gyrofrequency which is simply related to the 
Alfven velocity (3 a expressed in units of speed of light as 

Ba = A (15) 

a 

The phase and group velocities of these waves (in units of speed of light): j3 p h(kj) = uij/kj 
and /3 gr (kj) = duj/dkj, respectively, can be obtained from relations (|T2"D and (0). The 
parameter t p , which is a typical time scale in the turbulent plasma, is defined as 



^4«e(^T](?-l)^ (16) 



V 



2 e \B 2 /(87r 



The parameters important for our problem are Q e , a, q, k m i n and the ratio of plasma 
turbulent density to magnetic energy density; f tU rb = (87r£ tot / 'B 2 ). The above equations 
hold for both electrons and protons. In what follows we consider only interaction and 
acceleration of electrons. 



3.1. Critical Angles 

As described in DP, in general, four values of kj contributes to the Fokker-Planck 
coefficients, except for \x = and 7 < and for the pitch angles between the two critical 
values when only two values of kj are involved. For non-relativistic electrons one of these 
roots is due to resonant interaction with ion-cyclotron waves with kj ^> 1, so that its 
contribution to the coefficients is negligible {x{kj) <C 1). The main contribution to the 
Fokker-Planck coefficients then comes from the electron's interaction with whistler mode 
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and higher-frequency electromagnetic branch. The last interaction, allowed only in low 
a-plasmas (a < 1), occurs for pitch angles greater than some critical angle /i^™ at which 
the group velocity of electromagnetic wave becomes equal to the component of the electron 
velocity along the magnetic field, (3 gr = Similarly, there exists another critical angle 
\f^ r < fil™ for interaction of electrons with electron-cyclotron waves (the high-/c end of the 
branch which is commonly refered to as whistler branch at u{k) < 1), It is well known 
that for these angles the quasilinear approximation breaks down and the Fokker-Planck 
coefficients become infinite. This has a minor consequence for the general process considered 
here (Steinacker and Miller 1992, DP). However, in order to clarify some of the behaviors 
of the coefficients we give a brief description of these angles. 

The dependence of these two critical angles on electron velocity f3 and plasma 
parameter a are shown on Figure [TJ The upper and lower curves correspond to \f™ and 
respectively. Thus for < fi < y^. there are three roots coming from the interaction 
with forward (2 roots) and backward (1 root) moving electron-cyclotron waves. Similarly, 
for 1 > fi > if™ three are roots, two from the electromagnetic branch and one from 
the backward moving electron-cyclotron branch. But for \i e ™ > \x > only the later 
root exists. The fourth and unimportant root from the (backward) ion-cyclotron branch 
mentioned above is common to all three cases. Same situation holds for negative values of 
\i with the role of the forward and backward waves reversed. 

For each value of a there exist a critical velocity (5 cr or kinetic energy E cr below which 
\f™ = 1 and interaction of electrons with electromagnetic branch is not allowed. Figure |2] 
shows the variation of E cr with the plasma parameter a. For high values of a (high plasma 
density, low magnetic field) [3 cr — ► 1. In the opposite case of a — > ( very high magnetic 
field and/or low plasma density) the two curves merge to the solid line on Figure [I] which is 
described by the expression \i cr = (1 — y/1 — (3 2 )/ (3. 

The effect of this behavior can be seen on Figure |3] where we show the dependence of 
one of the Fokker-Planck coefficients, D pp , on electron pitch angle in a plasma with a = 0.2. 
As we approach one of the critical values of the pitch angle //^J 1 or //^ the coefficient 
becomes infinite as indicated by the sharp cusps. For each energy the first and second peaks 
occur at \x e ^ r and y^ 1 , respectively. 
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Fig. 1. — The dependence of the critical angles (lower curves) and (upper curves) 
on electron's velocity (3 for different plasma parameters a. For a = (the empty space) 
Vtr = VtT (solid line). 
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Fig. 2. — The plasma parameter dependence of the critical value of electron kinetic energy, 
below which interaction with the upper electromagnetic branch is not allowed. 

4. RATIOS R x AND R 2 . 

Using equations (Q) to ( |Tt]| ) the ratios defined in equation (^) become 

R = SjLi/^(fc 3 -) 2 x(fc 3 -) 

Y.UM k i)W-i*M k i))x( k i) (17) 
2 Ef=i(/?-^(%)) 2 x(%) ' 1 J 

which depend on the electron pitch angle and velocity, and on the plasma parameters a 
and q. Often, in particular for \f™ > \i > only one root has significant contribution. In 
this case the above expressions simplify to R\ — R\ — uj 2 /(j3k — ficu)^ 2 . 

Since the derivation of transport equation (U) requires these ratios to be much less than 
1 we want to determine the plasma conditions and energy range for which this requirement 
is satisfied. It is well known that relativistic electrons (and all protons), interact mainly 
with Alfven waves. In this case (see also §6 below) R\ oc (f3 a / ' (3) 2 <C 1 if (3 a 1 for 
relativistic electrons (or (3 f3 a for protons). This condition tends to breakdown as we 
go to the lower energies, where the interaction of the electrons with other plasma waves 
become more significant than the interaction with Alfven and ion-cyclotron waves. 
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Fig. 3. — The dependence of the Fokker-Planck coefficient D pp on pitch angle for electrons 
with two different values of kinetic energy in plasma with a = 0.2 and q = 1.6. For smaller 
values of a the flat, single root regime in the middle range of \i becomes smaller. This regime 
extends to the whole range of /i for larger values of a and non-relativistic energies. 



4.1. Electrons with 90° Pitch Angles. 

The above equations can be solved analyticly for \x = when the details of the 



distribution £(k) are not important. For this case the resonant condition (|P2|) is simplified 
to uj(kj) = ±1/7 and the ratio R\ = f3 2 h {k)/(3 2 . The ratio R2 = because D w = 0. 
Substitution of this in equation (|T3|) gives four symmetric roots (two from the Alfven branch 
and two from the whistler branch) which differ only in their signs of Uj or kj. We then have 
the analytic relation 

(1 + 'j)(Ya 2 + 7 - 1) 

In the extreme relativistic regime 7 ^> <5 _1 , R\ = -^r^ < 1 and for 1 <^ 7 <^ 
Ri = 7 _1 (a 2 + 1, but for non-relativistic electrons R\ ~ l/2a 2 which can exceed 

unity if a 2 < 1/2. Figure |] shows the relation between the electron kinetic energy and the 
plasma parameter for five different values of this ratio. As evident, for small a the condition 
(0) required for validity of equation (||) is violated up to very high energies (Note that for 
a « ~ 0.024 this is true at all energies). For example, for plasmas with a ^ 0.1 the 
diffusion approximation becomes valid only for electrons with kinetic energy exceeding few 
MeV. This means that for low energy electrons in low-a plasmas we cannot use equation 
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Fig. 4. — The dependence of the electron kinetic energy E on plasma parameter a for five 
different ratios R\ for electrons with pitch angle /x = 0. 

(f|) and we must revert back to the original Fokker-Planck equation ([I]). 

In the opposite case of R\ ^> 1 (low energies, a 2 < 1) the scattering time becomes 
much greater than the acceleration time so that the electron pitch angle will not change 
significantly in times of order of the acceleration time scale. We can, therefore, neglect the 
small and D m 



terms in equation (JJ) and obtain the simple diffusion equation 

df 



df 
dt 



]__d 

p2 Qp 



0. 



(19) 



Thus for this special case of /x = (i.e. 90° pitch angle) the Fokker-Planck equation reduces 
to the pure diffusion equation in momentum (or energy) space. ^From inspection of Figure 
f| we can see that this equation must be used up to several hundreds of keV for a < 0.3. 



4.2. Electrons with // ^ 

The situation is more complicated for electrons with /i ^ because the resonant 
condition ( p!2[ ) is no longer simple and there are multiple roots kj contributing in the 
summations in equations (|17D. Analytic approximation for the ratios is not possible in this 
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Fig. 5. — The dependence of the ratio R\ on electron energy for four different values of pitch 
angle in the plasma with parameter a = 0.1 and spectral index q = 1.6. 

case and we have to evaluate these ratios numerically. Figures [5| - § show variation of R\ 
with energy for different values of fi in plasmas with four different values of a. As evident, 
these ratios vary considerably with /i but in general they tend to increase with decreasing 
energy. Furthermore, they depend on the plasma parameter a and the spectral index q. 
The ratios R± and R2 tend to increase with decreasing a (high magnetic field, low plasma 
density) but they are not very sensitive to q (compare Figures |5]and|B|). The discontinuous 
behavior of the curves arise from the existence of the two critical pitch angles discussed in 
the §3.1. For given a and fi there exists a maximum of three and a minimum of one energy 
for which this value of fi is critical and leads to a discontinuity (except /i = 0, which is not 
critical for any energy). Unlike the Fokker-Planck coefficients, their ratios do not become 
infinite at critical values of /1 or energy. In §5 we will derive the asymptotic expression for 
the ratios R\ and R2 as a function of a for low energy electrons. 

Following the same argument used in the case of zero pitch angle electrons we can 
conclude that for low energy electrons the distribution function can be obtained from the 
solution of the following equation 




OF 



(20) 



dt 



dz 



pp dp 
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Fig. 6. — Same as Figure || with a = 0.1 and q = 3. 
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Fig. 7. — Same as Figure || with a = 0.3 and q = 1.6. 
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Fig. 8. — Same as Figure |5| with a = 0.8 and q = 1.6. 



where D£ p is given by equation (|10|) for each /i. The advective term (second on the left hand 
side) can also be neglected if the traverse time L/v is much larger than the acceleration 
time p 2 /D pp and if D pp and the source term are independent of z (homogeneous plasma). 



Equation (|20|) then becomes identical to equation (19) for zero-pitch angle electrons. Once 
the electron achieve higher energies and moves above the line R\ — 1 so that the conditions 
(0) are satisfied, then we revert to equation (^). In the transition region neither equation 
(||) nor (PDD are valid and the resulting equation for isotropic pitch angle distribution is 
almost as complicated as the original equation ([I]). 



NON-RELATIVISTIC APPROXIMATIONS 



For low energy electrons 7 — > 1 + /3 2 /2 and the resonance condition (|T2| ) combined with 
dispersion relation (O) gives 



fi(3k 3 - (- + 2fi 2 )(3 2 k 2 - fi(3{l - a 2 )k + a 2 = 0, 



(21) 
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where we have assumed a 2 3> j3. It can be shown then that 11^. = -^j^cn <C 1 and — > 1 
so that the single root solution is applicable for almost all pitch angles with the root 

kj ~ -/3~s a §/i-§, 1 > /j, > n™. (22) 



In the limit of \i — > equation fl2"2] ) has three roots one of which is very large (k — > oo) 
compare to the other two roots which are 



V2a 

T 



%-±^V> A*</£- (23) 



ec 



This expression is nearly identical to the value of kj obtained from equation (p2|) at ll — fi[ 
so that joining this two relations at /i = fi cr will give an approximate description for all ii. 

Now substitution of these values for kj in equation ( |T2"D gives the values of the resonant 
frequencies as 

u3 - 1 ~ \ + V%<V< 1, (24) 

Then to the first order the phase velocity of the wave is obtained from equations fl2~2"|) to 
([24]) to be = ^ ~ p. Combining this with the dispersion relation (|IB|) we obtain the 



group velocity 

^ = cL^2a2 f 2/i/5 + o(/3 2 ), ^ c </i<l, 
dfc A; 3 | /x < iL e ^ r . 

Using the above expressions we can now evaluate the Fokker-Planck coefficients and 
their ratios in the non-relativistic limit. For the single root case these ratios are simplified 
to Ri = R 2 , = uj 2 /{[3k — fiu) 2 and in the small range of ll < (i™ we can use the asymptotic 
relation of equation dl8|) . Substituting of equations (|22|) through (|25|) into equation ([17]) 
gives the non-relativistic values of these ratios as 

u" < a < 1 (26) 

These values agree well with the numerical results shown on Figures [5| to ||] at low 
energies. 
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Fig. 9. — Dependence of acceleration time r a = {p 2 / D pp ) in units of r p on pitch angle \i in 
plasma with a = 1 and q = 5/3 for electrons with different kinetic energies. 

5.1. Acceleration Time 

For non-relativistic electrons and especially for a < 1 the ratios R\ and R2 exceed 
unity and the acceleration process is described by equation (^). We can therefore define 
an acceleration time r a = p 2 /D pp . [Note that this is different from r ac = 1/k 3 defined for 
relativistic regime in combination with equation (Q); see below]. Figure || and show the 
dependence of this time (in units of r p ) on fi at three different values of kinetic energy in a 
plasma with q = 5/3 and a = 1 and a = 0.2, respectively. For a > 1 the critical pitch angle 
[i e ™ = 1 and the electrons interact mainly with the whistler and electron-cyclotron waves 
and not with the electromagnetic branch. Over the wide range of \i the resonant interaction 
occurs only with backward (for /i > 0) or forward (/i < 0) moving whistler waves. For 
a = 0.2 we have \f™ < 1 (Figure |T^) and there are three distinct regions on the plot with 
two discontinuous changes. Note that the acceleration time of electrons with pitch angles 
/1 > /i^ becomes large for all energies because of the dependence D pp oc (1 — fi 2 ). 

Using values of resonant wave vectors obtained in equation (||) we get the following 
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Fig. 10. — Same as Figure |9| with a = 0.2 and q = 5/3. 



expressions for the acceleration time 



Tail*) 



,(l-g)/3 2{q+2) 7-g 



l_. t a -Oi 3 P~ 



(27) 



To demonstrate the dependence of the acceleration time on electron energy alone we 

lera 

(r.) 



plot on Figure |TT| the averaged acceleration time 

2p 2 



(21 



From equations ( P7|) we also obtain an approximate analytic expression 



(T ) (2 + q) (8 + g) 2( 9 +2) 7-9 

r„ 18 



(29) 



where we have ignored the small contribution from electrons with \x > \f™ because in 
most cases \i e ™ is close to 1. This approximation agrees with 10% precision with numerical 



results shown on Figure |TT] for values of a > 0.3 and electron energies of order of keV. 



Simple analytic solutions are not possible for the plasma conditions when a significant 
fraction of electrons is able to interact with the electromagnetic branch. As can be seen 
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Fig. 11. — Dependence of averaged acceleration time defined for low- energetic electrons in 
units of r p (equation |28|) on kinetic energy of electron in plasmas with two different values 
of spectral index q and a = 1, 0.3 and 0.1 from top to bottom. 
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from Figure 11 in the case of small a the averaged acceleration time is a complicated 
function of energy and plasma parameters. In order to determine the highest energy for 
which the approximate analytic expression Q2"5|) is valid one should refer to Figure |2|. For 
example, for a = 0.3 we obtain E cr ~ 10 keV. This is in agreement with Figure |TT] which 
shows the correct power law dependence below 10 keV. 



6. RELATIVISTIC APPROXIMATIONS 

To complete our investigation of the Fokker-Planck equation (JJ) in application to 
the electrons in a turbulent plasma in this section we consider the case of the relativistic 
electrons for which Ri^(n,p) <C 1 and the transport equation (§) is applicable. For this 
regime the equations (|5]) and (|7|) define the average scattering time r sc and, now the longer 
acceleration time r a as 

Tsc = J_ ^(l-/! 2 ) 2 /^, (30) 

Tac = J^d^R.-RjY (31) 

Figures [T2| and Figure [13] show the variation with energy of these times (in units of 
t p ) obtained from numerical integration of the above equations for different values of the 
plasma parameter a and spectral index q of the plasma turbulence. The scattering time 
is similar to the scattering time r sc = 1/(Z) MM ) one would define in the non-relativistic 
regime. Thus the curves in Figure [12] provide a good estimate of scattering time scale at all 



energies. However, the acceleration time r ac in equation (|3T|) is different from the low-energy 



definition r ac = p 2 /D pp used in §5.1. These two definitions differ by the presence of R\ term 



in the denominator of equation ([31]). It turns out that for relativistic electrons in plasmas 
with small Alfven velocity R2 — R\ <C 1 so that R\ <C R\ and the two equations give very 



similar results. This can also be seen from comparison of Figures |TT] and |13| As we have 
shown in §4, this condition is not true in the non-relativistic regime so that equation ( |3"T| ) 
when extended to low energies is in error by several orders of magnitude. 



From figures 11, 12 and 13 we see several distinct regions of energy with different 



behavior of the curves. 

a) Extremely Relativistic Energies: 7 ^> S^ 1 (3 a hi (3~ l . 

The power-law energy dependence { r ) 2 ~ q ) of t sc and r ac seen in these figures at high energies 
is well known in the literature (e.g. Schlickeiser 1989; Schlickeiser et al. 1991). Relativistic 
electrons interact with Alfven waves which have a simple dispersion relation u = j3 a k. 



-19- 



100 



10 

Q. 

b 

e 1 



0.1 



0.01 



0.001 



0.0001 




0.001 



0.01 



0.1 



1 10 

E (MeV) 



100 



1000 10000 



Fig. 12. — Same as Figure |TT] except for r sc from equation ( [30]) and for a = 10, 1.0 and 0.1 
from top to bottom. 
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Fig. 13. — Same as Figure |TT] except for r ac from equation (|3ll) and for a = 10, 1.0 and 0.1 
from top to bottom. 
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Recall that 5 = m e /rrii and (3 a = \f8/a is the Alfven speed in units of speed of light, which 
we assume to be less than one. These electrons also interact with long- wavelength whistler 
waves which, to the first order in /3 a , have the same dispersion relation as Alfven waves. For 
both waves the resonant frequencies lor <C 5. 

Below we derive analytic expressions which show the above power-law dependence as 
well as the dependence of r sc and r ac on the plasma parameters a and q. Using the above 
simplified dispersion relation and the resonant condition ([TJ) we get two resonant values for 
the frequency and the wave vector: 

U 1 1 1 Pa „ x 

h ^ — r , ^ — =-, 32 

7 A* + Pa 1 H + Pa 

U 1 1 1 Pa ,„x 

k 2 ~ —, uj 2 ~ — . (33) 

7 A* - Pa II*- Pa, 

Note that the demarcation energy for this case comes from the condition ur 5 and is 
different from the usually assumed limit 7 3> 8 . Substitution of the above solutions in 
equations ( |30| ) and ( pT|) gives the asymptotic expressions for scattering and acceleration 
times 

T tA^ w ' Tr-^'wzMv (34) 

where 

Wa) = L \» + p a \* + \»-pr (35) 

and the power indices s sc — q — 1, s ac = 1 — g for scattering and acceleration times, 
respectively. 

One can show that in the case of small P a the function I s {P a ) simplifies to: 

f (2)9 a ) 1 -'/o 1 (1 - f) 'y 2 dt, *>1 
J s (/5 a )~ -l/21n(/5 a /4), s = l (36) 

i (T^fcy + ^a) 1 " 5 , 5<1 

Now we can use the above results in equations ( |34"D to obtain the dependence of relativistic 
times on plasma parameters. We will consider spectral indices in the range 1 < q < 4. 
In the case q > 2 the scattering time for electrons of relativistic energies has a power law 
dependence on plasma parameter, r sc oc a q ~ 2 . In the opposite case of q < 2 we have almost 



no dependence of scattering time on a. We can see these behavior in Figure for q = 5/3 
and q = 3. 

To the highest order in f3 a and for all q in the above range, r ac oc P~ 2 . Thus the 
relativistic acceleration time is proportional to a 2 and this dependence can be seen on 
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Figure [13]. The above results indicate that the ratio R\ ~ r sc / r ac has a strong dependence 
on a and q which leads to the fact that in low-a plasmas with spectral index q < 2 this 
ratio is of order of 1 even for relativistic electrons and equation is not valid. 

b) Intermediate Regime: 1 7 5 _1 /3 a ln/3" 1 . 

As can be seen from Figure [12] and O at 7 ~ 5~ l /3 a In f3~ l the curves begin to deviate 



from the described power-law because the interaction with whistler mode with the higher 
wave vectors begin to dominate. In this case 1 3> lor ^> 5 and the simple expressions 



such as (|32|) and ([34]) are not possible and one has to consider different combinations of 
plasma parameters in order to get analytic expressions for scattering and acceleration 
times correspondent to that particular plasma conditions. Different cases of interaction of 
electrons with whistlers have been considered in several papers (e.g. Steinacker and Miller 
1992, Hamilton and Petrosian 1992). 

c) Non-relativistic Energies: (7 - 1) < 1. 

Lastly, for non-relativistic electrons considered in §5 the resonant frequency ujr — > 1 and 
electrons interact most effectively with short- wavelength whistlers. As we described in §3.1, 
in plasmas with a < 1 electrons also interact with lower frequency electromagnetic branch 



of plasma waves. We see the corresponding change of behavior of the curves on Figures |TT 
and [12] as we approach non-relativistic energies. The interaction with electromagnetic waves 
is allowed only for electrons with energy greater than the critical energy E cr shown on 
Figure 0. In this regime the dependence of both scattering and acceleration times on energy 
deviates significantly from a power-law. For energies less than E cr equation (f28|) becomes 
valid and we again can see in Figure |TT] the power-law dependence of the acceleration time 
on kinetic energy of electron. 



7. SUMMARY 

The importance of the stochastic acceleration of high energy charged particles by 
turbulent plasma waves is well known. In this work we investigate the possibility of 
acceleration of the low energy background (often thermal) electrons by this process. We use 
the well known formalism developed over the years and specifically the formalism proposed 
by Schlickeiser and DP. In this paper we consider interaction of electrons with plasma 
waves propagating along the magnetic field lines. At all energies the plasma parameters 
which determine the acceleration rate of the electrons are the value of the magnetic field, 
the energy density of the plasma turbulence, the spectral index q of the waves and most 
importantly the plasma parameter a which is equal to the ratio of plasma frequency to 
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gyrofrequency of electrons. At low energies interaction with the electron cyclotron and 
electromagnetic branch of plasma waves become important. The second interaction is 
possible only when a < 1 and above the critical energies shown in Figure 0. The electrons 
which do not have enough energy to interact with the electromagnetic waves may still be in 
resonance with whistler and/or electron-cyclotron modes. 

From the above consideration we show that the ratio of the pitch angle to momentum 
(or energy) diffusion rates, or alternatively, the ratio of the acceleration to scattering times 
(equation [3D varies strongly with the plasma parameters and the pitch angle and energy 
of the electron. We give asymptotic analytic expression for these ratios showing their 
dependence on these parameters. In particular we show that this ratio becomes greater than 
one at lower energies which indicates that the usual transport equation derived for cosmic 
rays is not applicable for electrons for a wide range of energies. An approximate "critical" 
value of electron energy above which the dynamics of the electrons may be described by 
this transport equation is determined as a function of plasma parameter. We propose a new 
transport equation for non-relativistic electrons in low-a plasmas. 

We also show numerically and through asymptotic analytic expressions that when the 
above ratio is greater than one the acceleration time scale decreases with decreasing energy 
and could be very short for magnetic field of 100 G and turbulent energy density of less 
than 10 -4 of that of magnetic field energy density. To complete the discussion we also 
consider the diffusion rates for relativistic electrons for which the above mentioned ratio is 
less than one. We give analytic and numerical results on the acceleration and scattering 
time scales for different energies and plasma parameters. 

Based on these results we suggest the following scenario for acceleration of the 
background plasma electrons. The low energy (background) electrons can be accelerated by 
whistler, high-frequency electromagnetic or electron-cyclotron waves without a significant 
change in their pitch angle over time scale of order of the acceleration time r a = p 2 / D^/i). 
This process holds till electrons reach an energy where the scattering time becomes 
comparable to the acceleration time. Then intensive scattering along with acceleration takes 
place and for this stage we have to take into account all of the terms in the Fokker-Planck 
equation to describe the time evolution of the distribution function. When electrons 
accelerate to high enough energies (typically energies of tens MeV) their scattering time 
becomes much less than all other time scales, the distribution function becomes nearly 
isotropic and the well-known transport equation (M becomes applicable. For the two 
limiting cases with simple transport equations analytic expressions for scattering and 
acceleration times as a function of energy, plasma parameter and turbulence spectral index 
can be used. 
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In this paper we have considered interaction of electrons with plasma waves propagating 
along the magnetic field lines, in a "cold" background plasma. In subsequent papers we 
will consider waves with oblique propagation angles and include the effects of the finite 
temperature of the background plasma. 

We acknowledge J. Miller and P. Sturrock for useful discussions. This work was 
supported by NSF grant ATM 93-11888 and NASA grant NAGW 1976. 

REFERENCES 

Achatz, U., Droge, W., Schlickeiser, R., Wibberenz, G. 1993, J. Geophys. Res., 98, 13261 

Achterberg A. 1981, Astron.Astrophys., 97, 259. 

Beaujardiere, J.-F. de la and Zweibel E.G. 1989, ApJ, 336, 1059. 

Benz A. 1977, ApJ, 221, 270. 

Benz A. and Smith D. 1987, Solar Phys., 107, 299. 

Dung, R. and Petrosian, V. 1994, ApJ, 421, 550 

Hamilton, R.J. and Petrosian, V. 1992, ApJ, 398, 350 

Jaekel, U. and Schlickeiser, R. 1992, J. Phys., 18, 1089 

Jokipii, J.R. 1966, ApJ, 146, 480 

Kennel, C.F. and Engelmann, F. 1966, Phys. Fluids, 9, 2377 

Kennel, C.F. and Wong H.V. 1967, Plasma Phys, 1, 75 

Kirk, J.G, Schlickeiser, R. and Schneider, P. 1988, ApJ, 328, 269 

Lerche, I. 1968, Phys. Fluids, 11, 1720 

Melrose D.B. 1974, Solar Phys., 37, 353. 

Melrose D.B. 1980, Plasma Astrophysics, (Gordon and Breach, New York). 
Miller J. A. and Ramaty R. 1987, Solar Phys., 113, 195. 
Miller J.A. 1991, ApJ. 376, 342. 



-24- 



Schlickeiser, R. 1989, ApJ, 336, 243 

Steinacker, J. and Miller, J. A. 1992, ApJ., 393, 764 

Stix, T.H. 1992, Waves in Plasma, (AIP, New York) 

Sturrock, P.A. 1994, Plasma Physics, (Cambridge University Press) 

Tademaru, E. 1969, ApJ, 158, 959 



This preprint was prepared with the AAS IATeX macros v4.0. 



